Upgrading nirmatrelvir to inhibit SARS-CoV-2 Mpro via DeepFrag and free energy calculations.

The first oral drug for the treatment of COVID-19, Paxlovid, has been authorized; however, nirmatrelvir, a major component of the drug, is reported to be associated with some side effects. Moreover, the appearance of many novel variants raises concerns about drug resistance, and designing new potent inhibitors to prevent viral replication is thus urgent. In this context, using a hybrid approach combining machine learning (ML) and free energy simulations, 6 compounds obtained by modifying nirmatrelvir were proposed to bind strongly to SARS-CoV-2 Mpro. The structural modification of nirmatrelvir significantly enhances the electrostatic interaction free energy between the protein and ligand and slightly decreases the vdW term. However, the vdW term is the most important factor in controlling the ligand-binding affinity. In addition, the modified nirmatrelvir might be less toxic to the human body than the original inhibitor.

Machine learning combines atomistic simulations to predict SARS-CoV-2 Mpro inhibitors from natural compounds.

To date, the COVID-19 pandemic has still been infectious around the world, continuously causing social and economic damage on a global scale. One of the most important therapeutic targets for the treatment of COVID-19 is the main protease (Mpro) of SARS-CoV-2. In this study, we combined machine-learning (ML) model with atomistic simulations to computationally search for highly promising SARS-CoV-2 Mpro inhibitors from the representative natural compounds of the National Cancer Institute (NCI) Database. First, the trained ML model was used to scan the library quickly and reliably for possible Mpro inhibitors. The ML output was then confirmed using atomistic simulations integrating molecular docking and molecular dynamic simulations with the linear interaction energy scheme. The results turned out to show that there was evidently good agreement between ML and atomistic simulations. Ten substances were proposed to be able to inhibit SARS-CoV-2 Mpro. Seven of them have high-nanomolar affinity and are very potential inhibitors. The strategy has been proven to be reliable and appropriate for fast prediction of SARS-CoV-2 Mpro inhibitors, benefiting for new emerging SARS-CoV-2 variants in the future accordingly.

Correction: Rapid prediction of possible inhibitors for SARS-CoV-2 main protease using docking and FPL simulations.

[This corrects the article DOI: 10.1039/D0RA06212J.].

Correction: Characterizing the ligand-binding affinity toward SARS-CoV-2 Mpro via physics- and knowledge-based approaches.

Correction for 'Characterizing the ligand-binding affinity toward SARS-CoV-2 Mpro via physics- and knowledge-based approaches' by Son Tung Ngo et al., Phys. Chem. Chem. Phys., 2022, https://doi.org/10.1039/d2cp04476e.

Characterizing the ligand-binding affinity toward SARS-CoV-2 Mpro via physics- and knowledge-based approaches.

Computational approaches, including physics- and knowledge-based methods, have commonly been used to determine the ligand-binding affinity toward SARS-CoV-2 main protease (Mpro or 3CLpro). Strong binding ligands can thus be suggested as potential inhibitors for blocking the biological activity of the protease. In this context, this paper aims to provide a short review of computational approaches that have recently been applied in the search for inhibitor candidates of Mpro. In particular, molecular docking and molecular dynamics (MD) simulations are usually combined to predict the binding affinity of thousands of compounds. Quantitative structure-activity relationship (QSAR) is the least computationally demanding and therefore can be used for large chemical collections of ligands. However, its accuracy may not be high. Moreover, the quantum mechanics/molecular mechanics (QM/MM) method is most commonly used for covalently binding inhibitors, which also play an important role in inhibiting the activity of SARS-CoV-2. Furthermore, machine learning (ML) models can significantly increase the searching space of ligands with high accuracy for binding affinity prediction. Physical insights into the binding process can then be confirmed via physics-based calculations. Integration of ML models into computational chemistry provides many more benefits and can lead to new therapies sooner.

Searching for potential inhibitors of SARS-COV-2 main protease using supervised learning and perturbation calculations.

Inhibiting the biological activity of SARS-CoV-2 Mpro can prevent viral replication. In this context, a hybrid approach using knowledge- and physics-based methods was proposed to characterize potential inhibitors for SARS-CoV-2 Mpro. Initially, supervised machine learning (ML) models were trained to predict a ligand-binding affinity of ca. 2 million compounds with the correlation on a test set of R = 0.748 ± 0.044 . Atomistic simulations were then used to refine the outcome of the ML model. Using LIE/FEP calculations, nine compounds from the top 100 ML inhibitors were suggested to bind well to the protease with the domination of van der Waals interactions. Furthermore, the binding affinity of these compounds is also higher than that of nirmatrelvir, which was recently approved by the US FDA to treat COVID-19. In addition, the ligands altered the catalytic triad Cys145 - His41 - Asp187, possibly disturbing the biological activity of SARS-CoV-2.

Searching for AChE inhibitors from natural compounds by using machine learning and atomistic simulations.

Acetylcholinesterase (AChE) is one of the most important drug targets for Alzheimer's disease treatment. In this work, a combined approach involving machine-learning (ML) model and atomistic simulations was established to predict the ligand-binding affinity to AChE of the natural compounds from VIETHERB database. The trained ML model was first utilized to rapidly and accurately screen the natural compound database for potential AChE inhibitors. Atomistic simulations including molecular docking and steered-molecular dynamics simulations were then used to confirm the ML outcome. Good agreement between ML and atomistic simulations was observed. Twenty compounds were suggested to be able to inhibit AChE. Especially, four of them including geranylgeranyl diphosphate, 2-phosphoglyceric acid, and 2-carboxy-d-arabinitol 1-phosphate, and farnesyl diphosphate are highly potent inhibitors with sub-nanomolar affinities.

Insights into the binding and covalent inhibition mechanism of PF-07321332 to SARS-CoV-2 Mpro† † Electronic supplementary information (ESI) available: Includes a movie describing a representative binding process, the minimum distance between the non-hydrogen atoms of PF-07321332 and the non-hydrogen atoms of the catalytic dyad, the minimum distance dSγ–CN between the nitrile group of PF-07321332 and the sulfur atom of the Cys145 residue, the binding pathway of PF-07321332 to SARS-CoV-2 Mpro over the trajectories 0 and 20, the optimization starting from the ion pair Cys145−–His41H+, the optimized transition state for the hydrolysis of Int-3, and energy and Cartesian coordinates. See DOI: 10.1039/d1ra08752e

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been causing the COVID-19 pandemic, resulting in several million deaths being reported. Numerous investigations have been carried out to discover a compound that can inhibit the biological activity of the SARS-CoV-2 main protease, which is an enzyme related to the viral replication. Among these, PF-07321332 (Nirmatrelvir) is currently under clinical trials for COVID-19 therapy. Therefore, in this work, atomistic and electronic simulations were performed to unravel the binding and covalent inhibition mechanism of the compound to Mpro. Initially, 5 μs of steered-molecular dynamics simulations were carried out to evaluate the ligand-binding process to SARS-CoV-2 Mpro. The successfully generated bound state between the two molecules showed the important role of the PF-07321332 pyrrolidinyl group and the residues Glu166 and Gln189 in the ligand-binding process. Moreover, from the MD-refined structure, quantum mechanics/molecular mechanics (QM/MM) calculations were carried out to unravel the reaction mechanism for the formation of the thioimidate product from SARS-CoV-2 Mpro and the PF-07321332 inhibitor. We found that the catalytic triad Cys145–His41–Asp187 of SARS-CoV-2 Mpro plays an important role in the activation of the PF-07321332 covalent inhibitor, which renders the deprotonation of Cys145 and, thus, facilitates further reaction. Our results are definitely beneficial for a better understanding of the inhibition mechanism and designing new effective inhibitors for SARS-CoV-2 Mpro. The catalytic triad Cys145–His41–Asp187 of SARS-CoV-2 Mpro plays an important role in the activation of the PF-07321332 covalent inhibitor.

Insights into the binding and covalent inhibition mechanism of PF-07321332 to SARS-CoV-2 Mpro.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been causing the COVID-19 pandemic, resulting in several million deaths being reported. Numerous investigations have been carried out to discover a compound that can inhibit the biological activity of the SARS-CoV-2 main protease, which is an enzyme related to the viral replication. Among these, PF-07321332 (Nirmatrelvir) is currently under clinical trials for COVID-19 therapy. Therefore, in this work, atomistic and electronic simulations were performed to unravel the binding and covalent inhibition mechanism of the compound to Mpro. Initially, 5 µs of steered-molecular dynamics simulations were carried out to evaluate the ligand-binding process to SARS-CoV-2 Mpro. The successfully generated bound state between the two molecules showed the important role of the PF-07321332 pyrrolidinyl group and the residues Glu166 and Gln189 in the ligand-binding process. Moreover, from the MD-refined structure, quantum mechanics/molecular mechanics (QM/MM) calculations were carried out to unravel the reaction mechanism for the formation of the thioimidate product from SARS-CoV-2 Mpro and the PF-07321332 inhibitor. We found that the catalytic triad Cys145-His41-Asp187 of SARS-CoV-2 Mpro plays an important role in the activation of the PF-07321332 covalent inhibitor, which renders the deprotonation of Cys145 and, thus, facilitates further reaction. Our results are definitely beneficial for a better understanding of the inhibition mechanism and designing new effective inhibitors for SARS-CoV-2 Mpro.

Unbinding ligands from SARS-CoV-2 Mpro via umbrella sampling simulations.

The umbrella sampling (US) simulation is demonstrated to be an efficient approach for determining the unbinding pathway and binding affinity to the SARS-CoV-2 Mpro of small molecule inhibitors. The accuracy of US is in the same range as the linear interaction energy (LIE) and fast pulling of ligand (FPL) methods. In detail, the correlation coefficient between US and experiments does not differ from FPL and is slightly smaller than LIE. The root mean square error of US simulations is smaller than that of LIE. Moreover, US is better than FPL and poorer than LIE in classifying SARS-CoV-2 Mpro inhibitors owing to the reciever operating characteristic-area under the curve analysis. Furthermore, the US simulations also provide detailed insights on unbinding pathways of ligands from the binding cleft of SARS-CoV-2 Mpro. The residues Cys44, Thr45, Ser46, Leu141, Asn142, Gly143, Glu166, Leu167, Pro168, Ala191, Gln192 and Ala193 probably play an important role in the ligand dissociation. Therefore, substitutions at these points may change the mechanism of binding of inhibitors to SARS-CoV-2 Mpro.

Thermodynamics and kinetics in antibody resistance of the 501Y.V2 SARS-CoV-2 variant.

Understanding the thermodynamics and kinetics of the binding process of an antibody to the SARS-CoV-2 receptor-binding domain (RBD) of the spike protein is very important for the development of COVID-19 vaccines. In particular, it is essential to understand how the binding mechanism may change under the effects of RBD mutations. In this context, we have demonstrated that the South African variant (B1.351 or 501Y.V2) can resist the neutralizing antibody (NAb). Three substitutions in the RBD including K417N, E484K, and N501Y alter the free energy landscape, binding pose, binding free energy, binding kinetics, hydrogen bonding, nonbonded contacts, and unbinding pathway of RBD + NAb complexes. The low binding affinity of NAb to 501Y.V2 RBD confirms the antibody resistance of the South African variant. Moreover, the fragment of NAb + RBD can be used as an affordable model to investigate changes in the binding process between the mutated RBD and antibodies.

Benchmark of Popular Free Energy Approaches Revealing the Inhibitors Binding to SARS-CoV-2 Mpro.

The COVID-19 pandemic has killed millions of people worldwide since its outbreak in December 2019. The pandemic is caused by the SARS-CoV-2 virus whose main protease (Mpro) is a promising drug target since it plays a key role in viral proliferation and replication. Currently, developing an effective therapy is an urgent task, which requires accurately estimating the ligand-binding free energy to SARS-CoV-2 Mpro. However, it should be noted that the accuracy of a free energy method probably depends on the protein target. A highly accurate approach for some targets may fail to produce a reasonable correlation with the experiment when a novel enzyme is considered as a drug target. Therefore, in this context, the ligand-binding affinity to SARS-CoV-2 Mpro was calculated via various approaches. The molecular docking approach was manipulated using Autodock Vina (Vina) and Autodock4 (AD4) protocols to preliminarily investigate the ligand-binding affinity and pose to SARS-CoV-2 Mpro. The binding free energy was then refined using the fast pulling of ligand (FPL), linear interaction energy (LIE), molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA), and free energy perturbation (FEP) methods. The benchmark results indicated that for docking calculations, Vina is more accurate than AD4, and for free energy methods, FEP is the most accurate method, followed by LIE, FPL, and MM-PBSA (FEP > LIE ≈ FPL > MM-PBSA). Moreover, atomistic simulations revealed that the van der Waals interaction is the dominant factor. The residues Thr26, His41, Ser46, Asn142, Gly143, Cys145, His164, Glu166, and Gln189 are essential elements affecting the binding process. Our benchmark provides guidelines for further investigations using computational approaches.

Rapid prediction of possible inhibitors for SARS-CoV-2 main protease using docking and FPL simulations.

Originating for the first time in Wuhan, China, the outbreak of SARS-CoV-2 has caused a serious global health issue. An effective treatment for SARS-CoV-2 is still unavailable. Therefore, in this study, we have tried to predict a list of potential inhibitors for SARS-CoV-2 main protease (Mpro) using a combination of molecular docking and fast pulling of ligand (FPL) simulations. The approaches were initially validated over a set of eleven available inhibitors. Both Autodock Vina and FPL calculations produced consistent results with the experiments with correlation coefficients of R Dock = 0.72 ± 0.14 and R W = -0.76 ± 0.10, respectively. The combined approaches were then utilized to predict possible inhibitors that were selected from a ZINC15 sub-database for SARS-CoV-2 Mpro. Twenty compounds were suggested to be able to bind well to SARS-CoV-2 Mpro. Among them, five top-leads are periandrin V, penimocycline, cis-p-Coumaroylcorosolic acid, glycyrrhizin, and uralsaponin B. The obtained results could probably lead to enhance the COVID-19 therapy.